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PHOTOIONIZATION MODELS OF METAL-POOR EXTRAGALACTIC 
H II REGIONS AND THE PRIMORDIAL HELIUM ABUNDANCE 

V. Luridiana, 1 A. Pcimbert, 2 and M. Peimbert 2 
RESUMEN 

Discutimos el efecto de la contribution colisional a la intensidad de las lmeas de Balmer en la determination 
de la abundancia de helio primordial. Con este objetivo, presentamos un modelo de fotoionizacion de la region 
H II cxtragalactica de baja metalicidad SBS 0335-052. Mostramos que la abundancia de helio (Y) derivada 
en esta region depende de la contribution colisional a las lmeas de Balmer, tanto directa (a traves de una 
sub-estimation del cociente Hc/H real) como indircctamente (a traves de una sobreccstimacion de la extincion 
interestelar) . Dcsglosamos la forma en que ambos efectos afectan a la determination de Y. 

ABSTRACT 

We discuss the effects of collisional enhancement of Balmer lines on the determination of the primordial helium 
abundance. To this aim, we present a photoionization model of the metal-poor extragalactic H II region 
SBS 0335-052. We show that the derived helium abundance (Y) of this HII region depends on the amount of 
collisional excitation affecting the Balmer lines, both directly (through an underestimation of the actual He/H 
ratio) and indirectly (through an overestimation of the interstellar reddening). We detail how each of these 
effects affects the derived value of Y. 
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1. INTRODUCTION 

This contribution highlights the difficulties of ob- 
taining accurate determinations of the helium abun- 
dance (Y) in low-metallicity HII regions. High- 
quality Y values are necessary in order to reduce the 
uncertainty in the derived value of the primordial he- 
lium abundance (Yp ) , which is a fundame ntal quan- 
tity in cos mology: see also t he reviews bv iLuridianal 
(|2003|) and lPeimbertl l)2003l) for more complete dis- 
cussions of this topic and a quantitative estimate of 
the error budget in the determination of Yp . 

One of the main contributors to the overall un- 
certainty in Y for individual regions is the collisional 
enhancement of Balmer lines. In H II regions, the 
main process contributing to the intensity of Balmer 
lines is H + recombination, but a further contribution 
can arise from the collisional excitation of H° from 
the ground state. This contribution is generally a 
minor one, but it might reach a value of several per- 
cent under appropriate conditions; since it depends 
both on the fraction of H° and on the temperature, 
its evaluation requires both these quantities to be 
known accurately. 

To illustrate the difficulties inherent to this task 
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and suggest possible solutions, we present a pho- 
toionization model of SBS 0335-052, an extremely 
low-metallicity extragalactic H II region, and show 
how the abundance analysis of this object should 
be modified to take into account t he collisional en- 
hanc ement. In a forthcoming paper l|Luridiana et al.1 
l2003h we will present a more extended analysis 
of this object, as well as photoionization mod- 
els of two other low- metallicity regions, H 29 an d 
I Zw 18. See also [Davidson & Kinmanl lll985h 
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ISkillman fc Kennicuttl lll993l). IStasinska fc Izoto- 
(|200lf ). and lPeimbCTtTPeimbert.fcTuTidla^ 

for previous discussions of this topic. 

Low-metallicity regions are the most affected by 
collisions, since the collisional excitation of Balmer 
lines depends very steeply on temperature, which is 
especially high in these objects. As a rule of thumb, 
collisional contribution may be safely neglected in 
the analysis of objects with T e < 15, 000 K, while 
in hotter objects it may enhance the H/3 intensity 
by up to several percent. If not taken into account, 
this enhancement has two biasing effects on the de- 
termination of Yp: first, if the observed intensity 
is interpreted in terms of pure H + recombination, 
the He+/H+ abundance is underestimated; second, 
since the percentual collisional increase in Ha is al- 
ways higher than that in H/3, the measured Ha/H/3 
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ratio is larger than it would be without collisions, 
produ cing a spurious reddening in the Balmer spec- 
trum i Ferland fc Osterbrockl ll985V If this extra- 
reddening is mistaken for the effect of interstellar ex- 
tinction, the dercddencd ratios 7(A)//(H/3) derived 
for lines blueward of H/? are overestimated, and those 
derived for lines redward of H/3 are underestimated. 
The net effect of such bias is a decrease of the de- 
rived Y value, since the helium lines redward of H/3 
are globally brighter than those blueward of H/3 and 
weigh more in the abundance analysis. We will show 
this effect quantitatively in the case of SBS 0335- 
052, the second most metal-poor HII region known. 

2. SBS 0335-052 

SBS 0335-052 is a roundish region with a linear 
radius R ~ 80 pc. The reddening coefficient C(H/3) 
determined bv llzotov et alJ l|l999l) varies from 0.225 
to 0.33 along the slit; we will adopt in the following 
the value C(H/3) = 0.25, corresponding to the total 
F(Ra)/F(RP) in the slit. 

The model presented in this paper was c omputed 
with the photoionization code Cloudy 94 l|Ferlandl 
l2000albh : th e ionizing source was computed with 
Starburst99 IjLeitherer et We adopted 

the following input parameters: Q(H°) = 2.2 x 



10 53 s" 1 ; 



0.0010; t = 3.2 Myr; Z, 



gas 



0.0007; filling factor e = 1.9 x 10 3 ; covering fac- 
tor cf = 0.80; and a Gaussian density law N c — 
max(A^ o max exp(r/r )- 2 ,iV ( ! nin ), with A^ nax = 600 
cm" 3 , Ar™ in = 90 cm" 3 , and r = 225 pc. The 
results will be compared a gainst the observational 
data by llzotov et alJ l|l999|) . The model's output is 
presented both for the complete Stromgren sphere, 
and modified to take into account the a perture effect 
introd uced by the 1" x 5.4" slit used bv llzotov et alJ 
(|l999lh which was divided into 9 extractions (Fig- 
ure ^| . The model's predictions will be presented 
in the following figures as a function of the radius 
of the structure, both in linear (parsec) and angular 
(arcsec) units; the bottom axis is in arcsec for those 
plots that simulate observational (projected) quan- 
tities, and in parsecs for those plots that describe 
theoretical (radial) quantities. 

Figure [3 shows the profiles of the principal emis- 
sion lines predicted by our model. The theoretical 
points correspond to the simulation of the nine ex- 
tractions making up the sl it; they are compa red to 
the intensities observed bv llzotov et al 1 lll999h . 

Figure shows the radial (theoretical) behavior 
of various quantities related to the collisional contri- 
bution: from the top down, the ratio between the 
collisional and the recombination cmissivities of Ha 
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Fig. 1. Sketch of the slit used bv llzotov et alJ il999T> . 
and the nine extractions making it up, projected onto 
the model nebula (see text). The assumed distance is 
d = 57 Mpc. 
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Fig. 2. Predicted profile of some emission lines (solid line) 
compared to the observational data (dot-dashed line). 

and H/3; the differential contribution of the colli- 
sional and the recombination process to the Ha and 
H/3 luminosities; the cumulative collisional and re- 
combination Ha and H/3 luminosities; the hydrogen 
ionization structure; and, in the bottom panel, the 
electronic temperature and density. The top panel 
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Fig. 3. Predicted radial behavior of several quantities 
related to collisional enhancement: a) ratio between the 
collisional and the recombination emissivities of Ha and 
H/3; b) differential contribution of the collisional and the 
recombination process to the Ha and H/3 luminosities; 
c) cumulative collisional and recombination Ha and H/3 
luminosities; d) hydrogen ionization structure; e) elec- 
tronic temperature and density. 



of Figure El compares the projected profiles of the Ha 
and H/3 luminosities, with the aperture effect taken 
into account; the bottom panel of the same figure 
compares the predicted profile of the total Ha/H/3 
ratio to the profile of the Ha/H/3 ratio that would 
result from pure recombination. The recombination 
profile is almost constant, as expected from recom- 
bination theory for a structure with a fairly constant 
temperature profile; however, the collisional contri- 
bution, which peaks around R ~ 700 pc as a result 
of the interplay between the temperature and the H° 
fraction (Figure |3J), increases the total Ha/H/3 ratio 
along the whole diameter, and particularly in the 
outermost extraction (which is centered on R ~ 600 
pc). 

This increase in the total Ha/H/3 ratio would be 
observationally undistinguishable from an increase 
due to interstellar extinction: in Figure El we give 
examples of the reddening coefficient which would 
be observed along the nebula for various assumptions 
on the true extinction coefficient C(H/3) truo . 
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Fig. 4. Top: projected profiles of the Ha and H/3 lu- 
minosities, with the aperture effect taken into account; 
bottom: comparison between the predicted profile of the 
total Ha/H/3 ratio and the profile of the Ha/H/3 ratio 
that would results from pure recombination. 
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Fig. 5. Predicted reddening coefficient along the model 
nebula for various assumption on the true extinction co- 
efficient C(H/3) rea i. 

3. EFFECTS OF COLLISIONS ON THE 
DETERMINATION OF Y 

Table presents the values of some helium line 
fluxes, along with three sets of reddening-corrected 
line intensities. In the first set we neglect the col- 
lisional effects in the hydrogen lines; in the second 
we subtract the collisional contribution from the ob- 
served /(H/3), but we neglect the collisional effects on 
the observed reddening; finally, in the third we sub- 
tract the collisional contribution from both /(H/3) 
and C(H^) obs . The C(H/3) values associated to each 
intensity set are also listed. The table also presents 
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TABLE 1 

HELIUM LINE INTENSITIES AND 
ABUNDANCES 3, 



A 


F(\) 


1(A) 


/(A) 


/(A) 


F(H/3) 




WH/3) 


WH/3) 


4026 


1.294 


1.473 


1.490 


1.461 


4471 


3.404 


3.616 


3.659 


3.626 


4921 


0.821 


0.816 


0.825 


0.826 


5876 


11.476 


10.169 


10.290 


10.477 


6678 


3.237 


2.631 


2.662 


2.746 


Ha 


3.339 








C(H/3) 




0.25 


0.25 


0.21 


(y + ) 




0.0825 


0.0835 


0.0846 


AY 




0.0000 


0.0023 


0.0047 



a Observed He I line fluxes, relative to H/9, and corre- 
sponding dereddened intensities for three different as- 
sumptions on collisional effects (see Section |HJ. 



the resulting y + values, and the difference in the de- 
rived helium abundance AY for each of the three 
cases; these numbers are based on a very simple 
analysis, and are meant to provide only a rough in- 
dication of collisional effects. A more sophisticated 
analysis, based on better models and on an improved 
al gorithm for abundance analysis, will be presented 
in lLuridiana et all i|2003[K 

The following facts should be considered when 
reading this table: first, the effect depends on the 
number of lines used; second, the lines redward of 
H/3 are more brilliant, and therefore dominate the 
abundance determination; third, the effect might be 
a strong function of the position, especially in real 
nebulae, in which the behavior of the temperature 
and ionization degree are much less smooth than in 
a model. 

4. SUMMARY AND CONCLUSIONS 

In this work we discuss the influence of collisional 
enhancement of hydrogen lines on the accuracy of Yp 
determinations. The final g oal of this pro ject , whic h 
is presented in its entirety in lLuridiana et all (|2003tl . 
is to quantify the correction that should be applied to 
previous Yp determinations in which collisional en- 
hancement of Balmcr lines has been neglected, and 
to estimate the uncertainty attached to such correc- 
tion. In this framework, we also investigate to which 
extent the predicted collisions are model-dependent, 
and note that the occurrence of temperature fluctu- 
ations in real regions would enhance the collisions 
with respect to predicted values. 

We show that the helium abundance in individ- 



ual objects may have been underestimated as a con- 
sequence of the H/3 collisional enhancement. The 
underestimation in the Y values arises from the com- 
bination of two factors: the collisional H/3 enhance- 
ment and the collisional reddening. Since the last 
factor mimics a higher C(H/3), the quantitative ef- 
fect of the collisional bias on Yp is a function of the 
He I lines used in the analysis. 

As a consequence of this bias, the Yp value 
has also been underestimated. Since the effect 
on the individual Y values is more pronounced 
at high electronic temperatures, which characterize 
low-metallicity objects, the net effect on the Y — Z 
relation is a decrease of the slope. 
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